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In this paper, we study the physical meaning of the wavefnnction of the universe. With the con¬ 
tinuity equation derived from the Wheeler-DeWitt (WDW) equation in the minisuperspace model, 
we show that the quantity p(a) = for the universe is inversely proportional to the Hubble 

parameter of the universe. Thus, p(a) represents the probability density of the universe staying in 
the state a during its evolution, which we call the dynamical interpretation of the wavefnnction of 
the universe. We demonstrate that the dynamical interpretation can predict the evolution laws of 
the universe in the classical limit as those given by the Friedmann equation. Furthermore, we show 
that the value of the operator ordering factor p in the WDW equation can be determined to be 

p=-2. 

PACS numbers: 98.80.Qc, 98.80.Cq 


I. INTRODUCTION 

In quantum mechanics, the state of a particle is com¬ 
pletely described by its wavefnnction T'(a;, t). At the time 
t, the probability of finding the particle in an interval Ax 
about the position x is proportional to |lE'(a;, and 

thus p{x,t) = |T'(a:,t)/ is interpreted as the probability 
density and 'I'(a;, t) is called the probability amplitude for 
the particle. This is the statistical interpretation or the 
ensemble interpretation of the wavefnnction determined 
by the Schrodinger equation in the standard quantum 
mechanics. In quantum cosmology theory, the universe is 
described by a wavefnnction ip{hij, 4>) determined by the 
quantum gravity elation, called the Wheeler-DeWitt 
(WDW) equation [l|, 'H'4’{hij,(j)) = 0, where hij is the 
3d metric and / is a scalar field, rather than the clas¬ 
sical spacetime. In principle, the wavefnnction should 
contain all information about the universe 0, although 
it is hard to extract all the information from it At 
first glance, the wavefnnction of the universe should sat¬ 
isfy the statistical interpretation [^. But, if we adopt 
the statistical interpretation for the wavefnnction of the 
universe directly, it would be strange that the quantity 
/)/ denotes the probability density of finding a 
universe somewhere. 

The observed universe is unique. If we want to study 
the physical meaning of the wavefnnction of the universe 
with only one universe, the statistical or ensemble inter¬ 
pretation should be abandoned. In this paper, we study 
the properties of the wavefnnction of the universe in the 
minisuperspace model, in which the wavefunction of the 
universe can be determined by the unique parameter, the 
scale factor a to describe the time-dependent evolution of 
the universe. First of all, we find that the operator order¬ 
ing factor p representing the ambiguity in the ordering of 
noncommuting operators in the WDW equation should 
take value p = —2 due to the requirement of finiteness of 
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the wavefunction of the universe. Next, we show that the 
quantity p{a) = \ip(a)\‘^ for the universe is inversely pro¬ 
portional to the Hubble parameter of the universe, and 
represents the probability density of the universe staying 
in the state a during its evolution. We further show that 
the dynamical interpretation of the wavefunction of the 
universe can give inflation solutions of the small universe 
and the correct evolution laws of the universe in the clas¬ 
sical limit, as required by the correspondence principle. 
This paper is organized as follows. In Sec. n the WDW 
equation is applied to the minisuperspace model. Then, 
the formula of p{a) for the universe is obtained with a de¬ 
termined p. in Sec. EH The dynamical interpretation of 
the wavefunction of the universe is given in Sec. IIVI The 
correct evolution laws of the universe in the classical limit 
with dynamical interpretation is discussed in Sec. 0 In 
Sec. IVIl the formula of p{a) of the universe within the 
scalar field model is obtained. Finally, we discuss and 
conclude in Sec. IVHI 


II. WDW EQUATION IN THE 
MINISUPERSPACE MODEL 


Assumed to be homogeneous and isotropic, the uni¬ 
verse can be described by a minisuperspace model 
0 with one single parameter, the scale factor a. The 
Einstein-Hilbert action for the model can be written as 


S = 


/ Rc^ 
VlbTrG 


9£n / 

16c/ 




( 1 ) 


where the represents the energy density of the uni¬ 
verse, and the constant before £„ is chosen for later conve¬ 
nience. Since the universe is homogeneous and isotropic, 
the metric of the universe in the minisuperspace model 
is given by 

= cr^ . (2) 

Here, dQg = dr^/(l — kr^) + r^{d6‘^ + sin^ Odf/x^) is the 
metric on a unit three-sphere, N(t) is an arbitrary lapse 
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function, and cr^ = 2/37r is a normalizing factor cho¬ 
sen for later convenience. Note that r is dimensionless 
and the scale factor a{t) has length dimension From 
Eq. ([2]), one can get the scalar curvature 


7^ = 6 


6 


6k 


a'^c^N'^a a'^c^N'^a^ a'^a? ’ 


(3) 


where the dot denotes the derivative with respect to the 
time, t. Inserting Eqs. m cLiid- (|H|) into Eq. © , we can 
get 


f 6a‘^Nc^ / a^d ad'^ GenCi^ 

^ J IGttG ^ 


6a'^Nc^V f ( a}d aa 


IGttG J \N^c^ iV2c2 


_^ , GenO.^ 

+ + ka -^ 




dt, 


Nc^ 

In 


aa? , GenO? 




+ ka — 


dt. 


The Lagrangian of the bubble can thus be written as 


C = 


Nc^ /. ad? GenO? 


ka — 


2G V'"” c4 

and the momentum pa can be obtained as 


(4) 


dC 

Pa = ^ = -- 


c aa 


dd NG 
Taking N = 1, the Hamiltonian is found to be 

U = Pad - C 

1 f Gp^ c^ka o 

' “ + - ena^ 


2 V c^a G 


In quantum cosmology theory, the evolution of the nni- 
verse is completely determined by its quantum state that 
should satisfy the WDW equation. With = 0 and 
Pa = —a~^-^{a^-^), we get the WDW equation dii, 


-f )'(/'(a) = 0. (5) 


Here, k = 1,0,—I are for spatially closed, flat and open 
universe, respectively. The factor p represents the uncer¬ 
tainty in the choice of operator ordering. For simplicity, 
we take h=G = c= lui this paper. 


III. THE SQUARE MODULUS OF THE 
WAVEFUNCTION FOR THE UNIVERSE 

In order to study the physical meaning of the wave- 
function of the universe, we should discuss the square 
modulus of the wavefunction for the universe first. Since 
there is only one parameter a in Eq. ®, the complex 
function tpia) can be formally rewritten as 


■0(0) = i?(a) exp(i5'(a)), 


( 6 ) 


where R and S are real functions. Then the square mod¬ 
ulus of the wavefunction is 

p{a) = |0(a)|^ = R{af. 

It should be pointed out that, no matter what the phys¬ 
ical meaning it is, the value of p{a) should be finite. 

From E q. (IST). i t is easy to construct a conserved cur¬ 
rent as [III, [l^ 


Inserting Eq. 


f = -aP{ip*datp - 

daf = 0. 

^ into Eq. 0, we can get 
f = -oFR^S', 


(7) 

( 8 ) 


where the prime denotes the derivative with respect to 
a. Integrating Eq. ([8]) gives that 


aPR^S' = Co, 


(9) 


where the cq is a dimensionless integral constant. 

In the quantum Hamilton-Jacobi theory, the relation 
between the action and the momentum can be written 

i[ll, 


as 


dS dC 

Pa = = -aa. 

oa oa 


( 10 ) 


According to Eqs. and (HOD, we can get p{a) for the 


universe 


p(o) = R{aY = —- 


Co 


dP+^d' 

Using the definition of the Hubble parameter 


H{a) = -, 
a 


we can rewrite p{a) as 


p{a) = - 


Co 


aP+‘^H{a) ■ 


( 11 ) 


( 12 ) 


(13) 


In principle, the value of p{a) should be finite for any 
a, i.e., no matter whether the universe is very small or 
large enough. In both the inflation stage and the dark 
energy stage, the value of the Hubble parameter should 
be finite. From 


p{a —^ 0) = —- 


Co 


aP+^H{a) ’ 
we can get a boundary condition 

p + 2<0. 

On the other hand, the finiteness of p{a) with large a 


p{a —>■ oo) = — 


Co 


aP+2iJ(a) 
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requires that 


p + 2 > 0. 


The above two boundary conditions determine the value 
of the ordering factor to be p = —2. Thus, Eq. m is 
reduced to 


p{a) 


cq 

H{a)' 


(14) 


With the above formula at hand, we can discuss now the 
physical meaning of the wavefunction of the universe. 


IV. DYNAMICAL INTERPRETATION OF THE 
WAVEFUNCTION OF THE UNIVERSE 

In quantum mechanics, the wavefunction of 

a particle is interpreted as probability amplitude and 
|4'(a:, t)pAa; is the probability of finding the particle at 
time t in the interval Ax. This is the statistical inter¬ 
pretation or the ensemble interpretation of wavefunction. 
However, such an interpretation of wavefunction in quan¬ 
tum mechanics cannot be applied to the wavefunction 
of the universe in quantum cosmology. For an observer 
inside the universe, it is should be very strange if the 
wavefunction is related to the probability of finding a 
universe. So, the physical meaning of the wavefunction 
for the universe should be reinterpreted. 

From Eq. (HI, we can see that the value of p{a) for the 
universe only depends on the Hubble parameter H{a). 
p{a) is large when the universe expands slowly, and it is 
small when the universe expands quickly. So, it is ob¬ 
vious that p{a) for the universe represents the evolution 
velocity of the universe and thus relates to the dynamics 
of the universe. In this case, p{a) can be treated as the 
dynamical interpretation. Similar to the statistical in¬ 
terpretation of the wavefunction in quantum mechanics, 
the dynamical interpretation for the wavefunction of the 
universe can be explicitly described as: 

Jal is proportional to the time spent when 

the universe involves from the state ai to 02 . 

In this way, we have showed that the physical meaning of 
p{a) for the universe can be interpreted as the probability 
density of the universe staying in the state a during its 
evolution. 

Generally speaking, an interpretation for the physical 
meaning of wavefunction should satisfy the correspon¬ 
dence principle, i.e., the quantum cosmology can reduce 
to the classical cosmology in the classical limit within the 
dynamical interpretation. In fact, the dynamical inter¬ 
pretation completely depends on the evolution equation 
Eq. (IllL which we call the dynamical equation for the 
universe. As will be shown below, solutions of the WDW 
equation in the classical limit (a ^ 1) together with the 
dynamical equation can predict the same evolution laws 
of the universe as those given by the Friedmann equa¬ 
tion. The exponential expansion solutions of the early 


universe (a <C 1) can also be obtained from the WDW 
equation together with the dynamical equation m- 


V. THE EVOLUTION OF THE UNIVERSE 
WITH THE DYNAMICAL EQUATION 

If we want to probe the rationality of the dynamical 
interpretation of the wavefunction of the universe, we 
should verify whether the dynamical interpretation can 
give the correct evolution laws of the universe in the clas¬ 
sical limit or not [T^ . On the other hand, the quantum 
behaviors of the small enough universe should be reserved 
with the dynamical interpretation. The classical evolu¬ 
tion laws of the universe at different stages dominated by 
radiation, matter and dark matter, respectively, can be 
obtained by solving the Friedmann equation. The evolu¬ 
tion equations of the classical universe at different stages 
are shown in Tab. HI 


dominator 

density 

evolution 

radiation 

n=4 

a{t) <x (t + 

matter 

n=3 

a{t) oc {t + 

dark energy 

n=0 

a{t) oc. e* 


TABLE I: The classical evolution laws of the universe given 
by the Friedmann equation. The universe is dominated by 
different kinds of energies = Xn/aJ^ at different stages. 


Let us study the evolution laws of the universe with the 
dynamical interpretation of the universe. For simplicity, 
we only consider the case of the flat universe k = {). The 
WDW equation for the flat universe with energy density 
En can be written as 


Here e„ = Xn/d^, for n = 4,3,0 representing the uni¬ 
verse dominated by radiation, matter and dark energy, 
respectively la¬ 
in principle, the universe contains all kinds of energies 
at the same time, so the energy density should take the 
form of e = Eq + 63 + S4. In practice, the universe is 
dominated by one kind of energy e„ at a specific stage. 
During the evolution of the universe, n changes slowly 
from n = 4 to n = 0. For an arbitrary n, we can obtain 
the analytical solutions of Eq. (HU, 


V'n(a) 


104 Ju 


+C2W 


V 3-n/2 ) 

y^^3-n/2 y 

3-n/2 ) 


(16) 


Here, J^ixYs are Bessel functions of the first kind, 
Ya{xys are Bessel function of the second kind, and 
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First, consider the wavefunction of the WDW equation 
in the classical limit (a ^ 1). For x ^ — 1/4|, Bessel 

functions take the following asymptotic forms, 





VK12 — 7r/4), 


Yvix) 



VK12 — 7r/4). 


If the free parameters ci and ci in Eq. ra take the values 
of Cl = C 2 = c_ \J 7r/2, the wavefunction can be rewritten 
as 


n-2p-4 

Wn{ci) = c-a 4 exp 


3-n/2 


ie 


(17) 


where 9 = {3n — 2p — 16)7r/(4n — 24), and S' < 
0. From Eq. OT, we know that the wavefunction in 
Eq. (II3 describes an expanding universe as suggested 
by Vilenkin [HI- With the wavefunction above, we can 
get 


p{a) = |^/>(a)|" =c^a-P-2+"/2. 

Together with the dynamical equation (ED, the above 
equation gives 


d = 


-Co 

C^a-l+n/2’ 


(18) 


where cq < 0. Transforming the formula in Eq. ED into 

Cl 


and then integrating both sides of the equation, we have 


a(t) oc 




n 7 ^ 0, 
n = 0. 


The evolution laws of the universe from the WDW equa¬ 
tion in the classical limit (a ^ 1) are completely con¬ 
sistent with the solutions of the Eriedmann equation as 
shown in Tab. HI It is interesting that the evolution equa¬ 
tion of the universe derived from the WDW equation is 
independent of the operator ordering factor p, which def¬ 
initely means that p only represents the quantum effects 
of the universe (TB |. 

Next, consider the evolution of the early universe 
(a <C 1) within the dynamical interpretation. For x <C 
I— 1/4|, Bessel functions take the following asymptotic 
forms. 


J. 



1 

f(EE)’ 


y^(x) 





So the wavefunction Eq. ED for the small scale factor 
(a <C 1) can be rewritten as 


V'n(«) = c- 



^ i\n^^a^-P 

r(72)(6-n)'^\ 

V 2 \ 

. r( 72 -t l)(6-n)^ 

TrAn^^ / 


Here, we have assumed that p < 1, otherwise the wave- 
function is divergent. ci and C 2 take the same value 
c_ 1 / 71/2 as that determined by the solutions in the clas¬ 
sical limit. In this case, the probability density of the 
early universe is 


p(a < 1) = \ipnia < 1)1^ 


^2 r //(6 - n )^*^ 
27rA(( 


It is obvious that, when a ^ 1, p(a) approximates a 
constant denoted as po{n). When a <C 1, the dynamical 
equation ED can be rewritten as 


aP+^d = - 


Co 


Po{n)' 


which is solved by 


[-(p -f 2)co(t -f to)/po{n)] , 


p = -2, 


where H = —co/po{n). 

When the universe is very small, its behaviors are dom¬ 
inated by quantum effects and the evolution of the uni¬ 
verse depends on the operator ordering factor p: different 
p gives different evolution equations of the scale factor 
a{t). In fact, the ordering factor p has been determined 
by the boundary conditions of finite density matrix of 
the universe as p = —2. This specific value p = —2 
gives exponential expansion solutions for the early uni¬ 
verse, which is consistent with the result from quantum 
trajectory theory [3111. Thus, we conclude that the 
WDW equation together with the dynamical interpreta¬ 
tion can give exponential expansion solutions of the early 
universe and the correct evolution laws of the universe in 
the classical limit. 


VI. THE DYNAMICAL INTERPRETATION 
FOR THE MINISUPERSPACE MODEL WITH A 
SCALAR FIELD 


In quantum cosmology, the most widely used model is 
the minisuperspace model with a scalar field. For a FRW 
universe filled with a scalar field (f>, the WDW equation 
can be written as [3 [3 


'--aP- 
aP da da 


I 52 


U{a,<j)) 


^(a,</) = 0, 


where U{a,(j)) = o?{k — a?V{(j))). We can rewrite the 
wavefunction 7/(0, (/) as 7/>(a,(/) = i?(a, where 

both and 5(a,(/>) are real functions. Since the 


TT 















5 


wavefunction '0(a, 0) is a function of a and 6, the currents 
for WDW equation can be obtained as [l^ 


-aP{f*daip - fdaip*), 


-aPR^{a,p)daS{a,p), 

(19) 

- ipd^tp*). 


aF~^R^{a, f)d^S{a, p). 

(20) 


The quantum Hamilton-Jacobi theory gives the guidance 
relations Q 

daS{a, 4>) = -ad, 
d^S{a, (j)) = 

The currents and satisfy the continuity equation, 

daf + = 0. (21) 

Inserting Eqs. m cind- (j20l) into Eq. (EU, we can get 

daioFR^daS) - d^aP-^R^d^S) = 0. 

Applying the guidance relation to the equation above, we 
can obtain 

da{aPR^daS) + daiaP-'^R^d^S)a^<p/d = 0 . ( 22 ) 

Integrating Eq. (|2^ over a, we get 

aP+^R‘^d-aP+^R^<p‘^/d + J R^<pd{a‘^<p/d) =-cq. 

Let A{a, if) = R~^ f R^aP'^^ij){2(j)a + acj)aa)da, which gives 


P(a» = R^{a,(f) 


_^co_ 

aP+^a(l — a^ff) -\- A{a, cfj ’ 


(23) 


where = dip/da = cp/d, and paa = dpa/da. We can 
see when <p ^ 0, i.e., A(a, (/>) —^ 0 and —>■ 0, p{p, a) 

in Eq. (1^5)) will return to p{a) in Eq. (fTTl) . which indi¬ 
cates that the dynamical interpretation still holds for the 
wavefunction of the early universe with a slowly-rolling 
scalar field. 


VII. DISCUSSION AND CONCLUSION 


In summary, we have found the mathematical relation 
between the quantity p{a) for the universe and the Hub¬ 
ble parameter that p(a) is inversely proportional to the 
Hubble parameter H{a). We argue that p{a) is not the 
probability density of finding a universe somewhere, but 
represents the probability density of the universe staying 
in the state a. This presents a dynamical interpretation 
of the wavefunction of the universe. We have demon¬ 
strated that the dynamical interpretation can give the 
same evolution laws of the universe as those given by 
the Eriedmann equation, which satisfies the requirement 
of the correspondence principle that the quantum cos¬ 
mology should reduce to the classical cosmology in the 
classical limit. In this way, we have presented an inves¬ 
tigation of the physical meaning of the wavefunction of 
the universe. 

Another result is that the value of the operator order¬ 
ing factor p that represents the ambiguity in the ordering 
of noncommuting operators has also been determined. 
With the requirement of the finiteness of the wavefunc¬ 
tion of the universe, the ordering factor should take value 
of p = —2. This specific value of p can predict exponen¬ 
tial expansion solutions of the small universe. In fact, 
when the universe becomes large enough, the evolution 
of the universe is independent of the value of the order¬ 
ing factor p, which implies that p represents the rules of 
the quantization of the early universe and only dominates 
quantum behaviors of the universe. 


Acknowledgement 


Einical support from the National Key Basic Re¬ 
search Program of China (NKBRPC) under Grant No. 
2013CB922003 and National Natural Science Eoundation 
of China (NSEC) under Grant No. 61471356 is gratefully 
acknowledged. 


[1] B. S. DeWitt, Phys. Rev. 160, 1113 (1967). 

[2] J. B. Hartle and S. W. Hawking, Phys. Rev. D 28, 2960 
(1983). 

[3] A. Vilenkin, Phys. Rev. D 50, 2581 (1994). 

[4] S. Coleman, J. B. Hartle, T. Piran and S.Weinberg, 
Quantum Cosmology and Baby Universes, World Scien¬ 
tific Publishing Co. Pte. Ltd (1991). 

[5] N. Pinto-Neto and J. C. Fabris, Class. Quantum Grav. 
30, 143001 (2013). 

[6] N. Pinto-Neto, F. T. Falciano, R. Pereira, and E. S. San- 


tini, Phys. Rev. D 86, 063504 (2012). 

[7] S. P. Kim, Phys. Lett. A 236, 11 (1997). 

[8] S. Weinberg, Gravitation and Cosmology: Principles and 
Applications of the General Theory of Relativity, John 
Wiley and Sons, New York, 1972. 

[9] S. W. Hawking, Nucl. Phys. B 239, 185 (1986). 

[10] A. Vilenkin, Phys. Rev. D 33, 3560 (1986). 

[11] A. Vilenkin, Phys. Rev. D 37, 888 (1988). 

[12] L. P. Grishchuk, Class. Quantum Grav. 10, 2449 (1993). 

[13] A. Vilenkin, Phys. Rev. D 39, 1116 (1989). 




6 


[14] D. H. Coule, Class. Quantum Grav. 22, R125 (2005). 

[15] D. He, D. Gao and Q.-y. Cai, Phys. Rev. D 89, 083510 
(2014). 

[16] D. He, Q.-y. Cai, Science China 58, 079801 (2015); 


arXiv:1409.2115 (gr-qc). 

[17] D. L. Wiltshire, General Relativity and Gravitation 32, 
515(2000). 


